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Introduction and Background 



• Earth orbit and planetary (earth escape) missions necessitate precise control 
of propellants for long duration space environments 
- Cryogenic propellants necessitate use of accurate fluid management systems 


• Tank wall heat loads induce propellant warming and stratification 

• Vehicle maneuvers induce sloshing, increasing propellant warming and boil-off 

• Drives pressurization requirements and affects commodity mass 

- Current state of the art thermodynamic and thermal modeling typically performed 
independently 

- Models generally don't include fidelity to the level of 

• Boundary layer growth and energy/mass exchange 

• Stratification 

• Propellant slosh 

• Slosh baffles 

• Diffusion/evaporation 

• Gaseous helium infusion 

• Moving liquid vapor interface 

• Interactive thermal/thermodynamic coupling 




Model Description 



Fully Integrated models developed to produce high fidelity predictions for propellant 
thermodynamic states and characterization 


- Utilizes industry standard tools 

• Thermal/thermodynamic: Thermal Desktop and SINDA/FLUINT 

• CFD: FLUENT, FL0W3D 

- Can accommodate complex geometries, ullage pressurization systems, propellant 
conditioning, boil-off, boundary layer growth, slosh effects 

- Accounts for 

• Proper liquid/wall interface area 

• Proper liquid/vapor interface area 

• Development of "warm" layer or stratum 

• Proper mixing of fluid and ullage 

- Includes 

• Fluid conduction: stratification, boundary layer development 

• Convection: boundary layer development, propellant boiling 

• Mass transfer: diffusion, vaporization, condensation 

• Pressurization and venting 

• Dynamic liquid /vapor interface areas and liquid/wall interface areas 
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Flight and test data utilized to validate model predictions 
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Diffusion Modeling 
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• Initial conditions for T, P and gas mass fractions are set via calls to FLUINT subroutine (Changes 
thermodynamic states of the lumps ) 

• Mass flow rate connectors are used to manipulate amount of mixing adjacent control volumes 



Stratification Modeling 




• Stratification: "Temperature gradient 
within a fluid due to heat transfer by 
conduction and mass transport" 

For a container such as a tank, 
conditions at the wall are important 

- A temperature difference between 
wall and fluid 

- Heat transferred via conduction and 
fluid movement, 'boundary layer' 
(B.L.) 

- Fluid Behavior & fluid Properties 
Wall Material 

Modeling Requirements 

Sufficient resolution needed to capture 
stratification ( number of strata) 

Varying gravity/acceleration -> Local P &'h' 

Subroutines 

- Local boundary layer thickness 

- Mass flow rates 

- Heat T ransfer Coefficients 



Stratification Modeling: Simulation Results 



Predicted Boiling Heat Flux For Various Values of Wall Superheat And Gravity (Oxygen, 
Vertical Surface, Pressures = 15, & 30 psia, @ 0.0001 , 0.01 & 1 . q c ) 
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P = 15 psia, g/gc = 0.01 P = 30 psia, g/gc = 0.01 
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Wall Superheat (Twall -Tsat, Deg R) 


Percent of Bulk Fluid Thermal Stratification vs. Time For Various Values 
of Gravity (Liquid: Hydrogen. Oxygen and Nitrogen) 



Time (Minutes) 


• Full range of Convective heat transfer regimes accounted for 

• Boiling -> phase change and convective and radiation heat transfer 

• Boiling curves are unique for specific fluids 

• Boiling regimes -> phase change and convective and radiation heat transfer 

• Heat flux -> [T wa) | - T Satf | Uid ] and various wall surface effects & fluid 
properties 
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Stratification Modeling: Simulation Results 
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Cryo Propellant LOX Tank 
Tank fill level = 99% 
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• Higher g = more buoyancy forces, higher stratification (relative), 
eventually more mixing 

• Lower g = surface tension forces start having more of an effect, less of a 
buoyancy effect 






HELIUM 

SUPPLY 



• Sub cooling of a cryogenic propellant by 
helium injection is one of the most 
effective methods for suppressing bulk 
boiling 

• Implement finite rate heat transfer and 
instantaneous mass transfer model 

• Oxygen continues to diffuse into the 
helium bubbles until thermodynamic 
phase equilibrium is reached i.e. when the 
partial pressure of oxygen vapor in the 
bubble is equal to the saturated vapor 
pressure of liquid oxygen at the particular 
temperature of the liquid 

• The diffusion process (evaporation) 
results in a cooling effect caused by the 
vaporization heat of the surrounding liquid 
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Example Helium Infusion Results 



Helium Requirements for LOx Tank Pressurization 


a. 

< 



0.05 


0.1 

Helium Used [lbs] 


0.15 


29.5% Ullage 

32% Ullage 

35% Ullage 

48.5 Ullage 

63.5% Ullage 

— 66% Ullage 

— — 3% Ullage 
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Helium infusion is most effective with a small ullage percentage 
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Moving Liquid/Vapor Interface Modeling 



L/V Velocity vs. Flow Rate 
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L/V Velocity vs. Heiqth from Bottom 



Tank Height from Bottom (ft) 


Nominal Values Typical 
Cylindrical Section 

LH2 « 2.7 strata/min 
L02 « 1 strata/min 
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CFD to Thermal Mapping 
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•CFD fluid slosh 
simulation is carried out 
for the duration of the 
coast period using 
trajectory 6DOF data 
•Rotation and slosh 
events are identified 
and mapped from fine 
CFD volume mesh to 
coarse thermal surface 
mesh 

•Shown here is an 
example of a slosh wave 
in a LH2 tank during a 
coast period. 
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Droplet Modeling 



CFD Droplet Tracking 
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L 


• Receives a 3D VOF matrix from 
Flow3D as a simulation result 

• MATLAB subroutine parses the 
matrix. Uses an advancing front 
methodology to capture each 
droplet and their volumes and 
surface areas 

• Can be used as an input source 
for a slosh model 
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Droplet Modeling 



Droplet Evaporation 
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• Characterization of cryogenic 
droplet evaporation behavior 

• Result is a set of analytical 
equations to solve the amount of 
liquid evaporated during a slosh 
event for one droplet 

• Can be applied to multiple droplets. 

• Still needs to be verified 
experimentally, if possible. 
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Fully integrated models for generation of high fidelity predictions for 
propellant thermodynamic states and characterization have been developed 


- Models account for a variety of influences and reactions to the complexities of 

mission requirements in space environments 1 

- Models can be tailored to specific vehicles and mission timelines and durations 

- FLUINT pipe flow network simulation successfully used to model complex 
cryogenic thermodynamic behavior 

- Capability to include slosh baffle and propellant conditioning effects 


- Methodologies developed to map CFD predictions to thermodynamic networks 

within the models j 

- Option to include ullage and liquid droplet/vapor development 


Future Work 

- Dynamic simulation of liquid extraction during engine 

- Development of He usage due to dome wetting 

- Slosh induced ullage collapse modeling 



